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Abstract 



In the heavy quark hmit, the form factors for semileptonic and rare radiative B decays into Ught 
mesons are related by heavy quark spin symmetry. Here we compute the leading corrections of order 
A/mg to these symmetry relations, showing also how to include hard gluon effects systematically 
to any order in as{mQ). The subleading correction to the form factor relation for S — > vr can be 
computed exactly in the soft pion limit. For decays into vector mesons, the matrix elements of two 
local dimension-4 operators are needed, one of which vanishes in the constituent quark model. A 
few applications are briefly discussed. 
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A precise determination of the Cabibbo-Kobayashi-Maskawa matrix elements using ex- 
clusive heavy meson decays depends in a crucial way on our theoretical control over strong 
interaction effects in these processes. The heavy quark effective theory (HQET) provides 
a systematic framework within which the relevant matrix elements can be studied in an 
expansion in powers of I/tuq, the inverse heavy quark mass. We focus in this Letter on 
form factors for heavy-light transitions, which are relevant for semileptonic and rare decays 
such as -B ^ piui and B K*i~^i~. 

Heavy quark spin symmetry has been used to relate heavy-to-light form factors corre- 
sponding to different currents ||^, ^. These relations hold true in the kinematical region 
where the energy of the final light hadron E in the rest frame of the B meson is not too 
large. At large recoil E ^ ^qcd, a new symmetry comes into play, further simplifying the 
structure of the form factors [Q. At the same time, 0{as{mQ)) perturbative corrections due 
to the exchange of hard gluons become important. They can be computed in a systematic 
way in perturbation theory, and the corresponding 0{as) corrections have been evaluated 
inii- 

In both these kinematical regions there are corrections of order A/itlq to the symmetry 
relations for the form factors, arising from subleading operators in the HQET. The complete 
expansion of the individual form factors at order A/itlq involves in general a large number 
of both local and nonlocal contributions and appears to preserve little predictive power. 
Such analyses have been in given in for B ^ P and in ^ for the B ^ V form factors, 
where P and V are pseudoscalar and vector light mesons, respectively. We show in this 
Letter that, with certain modifications required to include hard gluon corrections, the 
symmetry relations can be simply extended to subleading order. The subleading corrections 
have a simple form and can be expressed only in terms of lowest order form factors and the 
matrix elements of two dimension-4 local operators. 

The hadronic matrix elements relevant for B decays into a light meson are parametrized 
in terms of form factors. There are three form factors for decays into a pseudoscalar, relevant 
for semileptonic decays (e.g. B — > vrez/) 

{P{p')\q^,b\B{p)) = U{p, + p',) + f4P,-p',) (1) 
and rare decays (e.g. B — > Ke^e^) 

{P{p')\qa,,b\B{p)) = is[{p^+p'^){p,-pl) - (p^+pI)(p^-p'^)] (2) 

The B decays into a light vector meson are parametrized by a total of seven form factors, 
four of which appear in semileptonic decay matrix elements 

{V{p',e)\q^^b\B{p)) = ig{q^)e^.x.el{p + p')x{p-p')a (3) 
{V{p\e)\q^,^,h\B{p)) = fiq^; (4) 

+ a+ (g2) (e* .p)(p + p')^ + {q^){t* ■ p){p - p')^ . 

and three form factors in rare radiative decays 

{V{p',e)\qa^^b\B{p)) = g+{q^)e^^xa^l{j> + p')„ + g^{q^)e^^xa^l{jp - p')„ (5) 

+Kq^)e^^xa{p + p')\{p - p')a{e* ■ p) . 
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While the vector and axial form factors are renormalization scale invariant, the form factors 
of the tensor current have a nontrivial scale dependence. Throughout in the following it will 
be understood implicitly that the scale used to define these form factors is /x = nib. 

Counting powers of ruf, coming from kinematical factors and the usual relativistic nor- 
malization of the \B{p)) states gives the heavy quark mass scaling properties of these form 
factors 

f+{y) + f-{y)'^m^^^^, -/-(y) oc s(?/)ocm^^/^ (6) 

and 

f{y) ccml^'^ , g{y) cc m^^^'^ , a+{y) - a^{y) cc m^^^'^ , a+{y) + a-{y) cc m^^''^ (7) 

-^-(y) oc mp, g+{y) + g_{y) (xn\^''^ , /i oc m^^/^ (8) 

These scaling relations hold at a fixed value of the light meson energy in the rest frame of B, 
Ey = mvy- Therefore for a consistent rrib power counting, we will use y = p-p'/ {msmv) as 
argument for the form factors, instead of q^. These parameters are related as = m\+my — 

2mBmvy- The range of y probed in semileptonic D decays is 1 < y < y^'^'^ = ^m^mv ~ ^■'^^ 
whereas in B decays the corresponding range is 1 < y < 3.0. 

For values oi y = v ■ v' not too far away from the zero recoil point ?/ = 1, heavy quark 
spin symmetry can be used to relate some of the form factors in (|l|)-® P, There is one 
relation among the B ^ P form factors 

(P-1) : U{y) - f4y) - 2mBs{y) = 0{m-'^') (9) 
and three relations among the B ^ V form factors 

(V-1) : - 9-iy) = -'^^BQiy) + c^K"'^') (10) 

(V-2) : g^{y) + g_{y) = 2myygiy) + — /(y) + Oim-'^') (11) 

rriB 

(V-3) : a+{y) - a^y) = 2g{y) - 2mBh + 0{m^'^') . (12) 

In the following we will review the derivation of these relations using a method which will 
allow us to include hard gluon corrections to all orders in as{mb), and to compute the 
leading 0{A/mb) corrections to them. 

We illustrate the principle of the method on the example of Eq. (^, and start by consid- 
ering the two currents = q'y^b and = iqcr^^v'^b. Both of them are matched in HQET 
onto operators of the general form 

Ji') = g-rWfe = cW {f^)q^^K + 4^ mv,K + ^ E (Z^)^; , (13) 

^ i=l 

with the most general set of dimension-4 operators with the same transformation prop- 
erties as the corresponding QCD current. Restricting ourselves to operators which do not 
vanish by the equations of motion, a complete basis can be chosen as 

Jl = q-if,ipK , Jl = qv^^iphy , J"^^ = qiD^K , (14) 

= V ■ id{q'Ji,h^) , Jl = v- id{qv^,K) , = id^{qh^) . 



3 



The Wilson coefficients of the dimension-3 operators Co\(m5) at the matching scale Ci{mb) 
are known to NNL order |T^. For they are given by 

co(mb) = 1 - - (11.40 - 0.79ne + 0.09) ( —) ' (15) 

ciim) = + (7.30 - 0.35n, + 0.02) f— V , (16) 



and for J' by 



c[)(mb) = -c;(mb) = -(1 - - (16.09 - 1.13n^ + 0.13) (—) ) . (17) 

The corresponding Wilson coefficients of the dimension-4 operators are given at tree level 
by 

6i(m,) = l, 62-6(^6) = (18) 
b[{mb) = b'^inib) = 1 , &3_6("^b) = . 

In addition to the local corrections coming from the matrix elements of the operators, 
there are also nonlocal subleading l/mj, operators which appear as time-ordered products 
of the leading dimension-3 currents with the l/m^ terms in the HQET Lagrangian 

-^^kin = K{iDfh^ , £mag = ^KgcTf^^F^^hy . (19) 

Reparametrization invariance |^ gives the relation Co(/i) = &i(/i) among the Wilson 
coefficients of the current J^, which should be satisfied to all orders in a^. No such constraints 
exist among the Wilson coefficients of J'^, for which reparametrization invariance has been 
broken explicitly by defining it in terms of t;^ (the corresponding constraints on the Wilson 
coefficients of cafj^^b are not sufficiently predictive to relate b'^ and c-). However, the condition 
v^J'^ = gives the following constraints for c^(yu) and &-(/i) separately 

cM = -c[{fi), b[{fi) = b',{fi), b',{fi) + bM + b',{fi) = 0. (20) 

The symmetry relation (^ follows from taking an appropriate linear combination of 
matrix elements of and between hadronic states (7r(p')| + ki J^|_B(p)), with ki to be 
chosen as explained below. In the full theory this matrix element can be expressed in terms 
of the form factors (|l]), (^ 

{nip')\J, + K,J'^\Bip))= (21) 
[/+ - Ki{mB - V ■p')s]{p + p')i, + [/_ + Ki{mB - V ■ p')s\{p-p')^ 

The same matrix element is given in the effective theory by 

{■k{p')\J^ + K.iJ'^\B{p)) = (co + fi:iCo)(g7^/i^) + (ci + Kic[){qv^K) (22) 



+^ \ H Ca(^)(co(/i) + /tiCo(^))(r{g7^/i„, iCa}) 

6 

+Ca(/i)(ci(/i) + Kic[{iJ,)){T{qv^hy, iCa}) + + '^1^0(^/1) 



i=l 
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At tree level, the combination of Wilson coefficients cq + kiCq multiplying the matrix 
element of q'^^h^ vanishes (with ki = 1). Since the contribution of qv^h^ is proportional to 
p^, it follows that the combination of form factors multiplying p'^ in ( pT] ) is suppressed by 
K/rrii,. This was the original observation of leading to the symmetry relation (P-1) 
From (|15D and (p^ one can see that cq + Cq = is accidentally satisfied also at one-loop 
order, and a correction appears ffist at two-loop order. We will take in the following 

= = 1 + (4.69 - OMn, + 0.04) ( (23) 

such that this cancellation persists to any order in as{mQ). Furthermore, this ensures 
also the cancellation of the nonlocal power-suppressed terms {Ca,qifj.hy)- The remaining 
dimension-4 local terms can be easily expressed in terms of leading order form factors plus 
the additional form factors 6±{y) defined as 

{7rip',e)\qtD,K\Biv)) = 5+{y){p, + p'^) + 5^{y){p, - p'^) . (24) 

The equation of motion for the field gives one constraint among these form factors {mB + 
v-p')S+{y) + {mB — v-p')S_{y) = 0, such that only one of them is independent. Inserting this 
into (|2^) and comparing the coefficient of on both sides, one finds the improved version 
of the symmetry relation (|]) 

(P-1') : /_(y) + K^mssiy) = - — (5+ - 5-) - (A - k,v ■ p')s{y) + 0{mf/') (25) 

rriB 



We used on the right-hand side of (p^) the tree level values of the Wilson coefficients 
bl'\mi,) (p!5D, i^n\j, which is sufficient for the accuracy intended. The corrections to (^) are 
suppressed by ^ relative to the leading terms in (|]), and come from matrix elements of 

dimension-5 operators which were neglected in 



Similar improved relations can be proved for the B ^ V form factors. The ffist symmetry 
relation (0) follows again from taking the matrix element 

{p{p',^)\Jf^ + i^iJ'f.lBiv)) = -(2g + Ki — (^+ -^-)) ie^,u\a^lpxp'„. (26) 

The same quantity can be expressed in the effective theory with the help of (|13D as 

{p\J^ + K^J'^\B) = (co(/i) + /tic[,(/i))(g7^/i,) (27) 



+ ^ \ H Ca(/i)(Co(^) + KiCo(/i))(T{g7^/l„, iCa}) + + 
^ ' a=kin,mag i=l 



The matrix element of qv^^h^ vanishes for this case by Lorentz invariance. The nonlocal 
l/rrih contributions cancel as before, leaving only local dimension-4 operators. They can be 
easily expressed in terms of leading order form factors plus one additional form factor 'D{y) 
defined as 

{V{p', e)\qtD,b\B{v)) = V{y)te,,xaelpxP^ . (28) 
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From ( ]27|) one finds now the following result for the corrected symmetry relation ([T0| ) 

^^i{9+ - 9-) + 2mB(7(y) = - Ug+ + g^) + —{g+ - g^) + 2I?(|/) ) + Oiml""^) . (29) 

Using the leading order symmetry relations (|ll]), (|12D on the right-hand side, this can be 
put in a form similar to (10) 

(V-1') : K,{g+ - g^) + 2mBg = ~2{mvy - A)g{y) - —f{y) - 2V{y) + O(m-=^/')(30) 

niB 

The remaining two relations ([TT|) and (0) can be derived in a similar way, starting with 
the two currents Js^ = qicr ^uV'^ l^b and JL = {g^u — v^v^)q'y'^'y5b. Both these currents can 
be matched onto HQET operators as in ([13D 

1 ^ ~ 

^ i=l 

where the dimension-4 operators appearing on the right-hand side are defined analogously 
to (|1^ with the substitution hy 75/^1,. The constraints v^J^^ = v'^J'^^ = give relations 
among the Wilson coefficients analogous to (|20|) 



4'^ M = if^) , + b^^ = , W if^) + W (/^) + ^ ^ (/^) = (32) 

The matching conditions for the current J^^ = qicr^uV'^lbb are 

co{mb) = di{mb) = 1 - - (16.09 - l.lSn^ + 0.13) (—Y (33) 

bi{mb) = 1 Mimb) = -1 , b3-Q{mb)=0, 
and for J'^^ = {g^^ - v^Vy)qYl-,b 

Co(mfe) = g;(mfe) = 1 - (10.88 - 0.77n^ + 0.11) (34) 
6^(mb) = b^inib) = -1 , b's^eimb) = . 

;;{') 



The values of Cq ^(mfe) quoted above correspond to the 't Hooft-Veltman scheme used in |]T3 



The two symmetry relations (|TT]) and (|T^) are obtained by taking again a linear combi- 
nation of the matrix elements of Js^ and Jg^, with to be determined below 

{p{p', e)| Js^ - K^J'^^\B{p)) = {g+{mB + v ■ p') + g-{mB - v ■ p') - n^f) (35) 
+ {v ■ e*){pf, + p'^) {-g+ - mBh{mB - v ■ p) + ^^f^5 [f - mB{mB - v ■ p'){a+ - a_)] 

+ {v ■ e*){p^ - p'^) [^9- + mBh{mB + v ■ p') + ^^^^5 [/ + rriBimB + v ■ p'){a+ - a_)] 

This relation involves a subtle point which deserves special attention. Usually, the matrix 
element (p|J5^|-B) is computed in terms of the form factors of the current ga^z^fe with the 
help of the relation cr^v75 = f e^z/a/jcr"^ (corresponding to the convention = 1). Beyond 
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tree-level, this relation is only true in the 't Hooft-Veltman scheme for 75. In any other 
scheme, there are calculable 0{as) corrections to the relation between the matrix elements 
(QCTf^i^b) and {qCfiulnb) , which cancel the scheme dependence in the individual Wilson coeffi- 
cients Co'\(mfe). For this reason, the use of the 't Hooft-Veltman scheme will be understood 
throughout in the following. 

The combination of matrix elements (^) is given in the effective theory by 

{p{p',e)\J5f, - KgJg^li?) = (co - K5C0) {{qif,l5hv) + {qVf,l5hv)) (36) 



' a=kin,mag i=l 



Using the tree level values of the Wilson coefficients, the contribution of the dimension-3 
operators vanishes with ^5 = 1, which gives the leading order symmetry relations (|lT|) and 
(|1^). The same cancellation will hold true to all orders in as{mQ) provided that K5 is chosen 
as K5 = Co(mfe)/c[)(mft) = 1 - as{mi,)CF/7!- + 0{al). 

The remaining matrix elements of dimension-4 operators can be expressed in terms of 
leading order form factors (0), (|^) and 3 additional form factors T>i,T>± defined as 

{Vip', e)\qiD,^,K\B{v)) = V,e; + V^{e* ■ p){p, + p'^) + V^e* ■ p){p, - p'^) . (37) 

The equation of motion for the heavy quark field iv ■ Dh^ = implies a relation among these 
form factors, such that only two of them are independent. 

Using standard HQET methods, one finds the following two generalized symmetry rela- 
tions, including subleading corrections 

(9+ + 9-)mB + v ■ p'{g+ - g^) - K^f = (38) 

{g+ + g^){A-v p') + —{g+ - g^){kv ■ p' - ml) + —V, + 0{mf/') 

niB ruB 



and 



g+ — g^ + 2m^h + K,5mB{a+ — aJ) = (39) 

— ig+ - g-) - 2mBh{K - v ■ p') - 2{V+ - D_) + 0(m,"'/') . 
rriB 

The ffist relation can be put in a form similar to (11) by inserting the lowest order result 
(0) for g^ — g_. One obtains in this way the improved version of the symmetry relation 



n]), including subleading K/rrib and hard gluon corrections 



(V-2'): g+ + g--2mvyg-n^—f = (40) 

TUB 

-2^(Ay - my)g{y) + — / + \{mvmByV{y) + Vr{y)) + 0(m,-'/') 
niB ttlb ttlb 



The analog of the symmetry relation ([12|) is obtained in a similar fashion from (|39D and 
reads 

(V-3') : K,{a+ - a.) + 2mBh - 2g = 2{^y - 2—)g{y) + -\f{y) (41) 

TflB TTIb ^B 

-2h{A - mvy) + ^my) - T^+{y) + T^-{y)) + (^«'^') 
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In conclusion, the leading A/mj, corrections to the symmetry relation (|^) ioi 13 ^ P 
decays are expressed in terms of one additional form factor S+ — S-, while the symmetry 
relations ioi B ^ V decays (|10D-(|12|) require three additional form factors V^y), Vi^y) and 
T)j^{y) — V^{y). All these four quantities are matrix elements of the two local operators in 
the effective theory qiD^{^^)hy. The relative simplicity of this result can be appreciated by 
noting that 22 new matrix elements, both local and nonlocal, are required for a complete 
1 / mq expansion of the individual form factors m B decays |^ . We note also that the 
local nature of these corrections should make their computation on the lattice feasible. In 
the following, we will estimate them using various approximation schemes and models. 



The form factors 5±(?/) (|2J) relevant for the i? — > P decays can be computed in the 
soft pion limit if P is one of the members of the Goldstone bosons octet. We use for this 
purpose the chiral perturbation theory for heavy hadrons developed in [Q. To lowest order 



in the chiral expansion, there is a unique operator which is a realization of qViDfJi^. Such 
an operator must satisfy the following conditions: a) transforms as (Sj,, 1r) + (Ii,, ^r) under 
SUl(?>) X SUr{?>), b) transforms in the same way under heavy quark spin rotations as qVh^ 
and c) vanishes upon contraction with v^. The only operator with these properties is 

qViD^K ^ ^«Tr [(7^ + v^)VH,{il + ^^a)] + [{^, + v^)-i,VH,{il - i.^)] , (42) 

where the Goldstone bosons are described by the field ^ = exp(j-M), and the heavy mesons 

B^^ are contained in the superfield Ha defined as in The low-energy constant a can 
be determined by taking the 5-to- vacuum matrix element of this operator with F = 7'^75, 
which gives a = — |/b^bA. 

The form factors S±{y) introduced in ( p4D can be computed now to the leading order in 
the chiral expansion. There are two contributions involving the operator (^21), the direct and 
the pole graphs. Their computation gives for the 5^ 7r+ form factors (this is multiplied 
with l/y/2 for the B~ 77° case) 

_gAfB TUB _ (0.5±0.1) GeV^ 

vp' + A- vp' + A ^^^^ 

with g the BB*it coupling and A = m^* — the hyperfine splitting in the B meson system. 
This matrix element can be also extracted from the current algebra calculation of 0. We 



used in the numerical estimate the recent CLEO measurement [T^ g = 0.59 ±0.12, A = 350 
MeV and /s = 180 MeV. 

For the B ^ V case, no predictions can be obtained from chiral symmetry. Some 
information can be obtained however from the constituent quark model. We will show in 
the following that the form factor T>{y) defined in (^) vanishes exactly in the quark model, 
which suggests that its real value might be very small. To see this, note that in the quark 
model the vector meson state is in general a mixture of S-wave and a D-wave states (^5*1 and 
^Di in spectroscopic notation ^'^"'"^Lj), whereas the B meson state contains only a 5- wave 
component ^Sq 

|V^(O,m, = +l)) = C505(|^|)|TT) (44) 
+CDM\k\) (^-^Y^MkM ii) + ;7^^2,+i(fc) ■ ^1 Ti + it) - 7f^2,o(^)| TT)) 

\B{v)) = M\k\)^\n-U) (45) 
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'D{y) is given by the product of an orbital overlap integral and a spin matrix element. The 
latter is only nonvanishing for the D-wave component of the vector meson. However the 
total result vanishes because of the orthogonality of the spin wavefunctions 

V{y) (X {V{q,m, = +l)|(e+i x q) ■ k\B{v)) (46) 



oc 



d'kM\k + q\)Y2:+iCk){e+i ■ k)M\k\) x (Ti + it | Ti - it) = . 



At large recoil Ey ^ ^qcd, the form factor T>{y) can be also computed in a light-cone 
expansion in 1/Ey. This kinematical region is at best appropriate io B ^ V decays at 
low values of g^, which is strictly speaking outside the validity domain of (pll). Still, such a 
calculation is instructive as an order of magnitude estimate. In this approach, the matrix 
element ( |25| ) is given by a diagram with a hard gluon exchanged between the current and 
the spectator quark. We will consider in the following the B ^ p case, for which one obtains 
at lowest order in 1/Ep 

^^^''^ = E^f'f^ Jo ^^^^^ I ^"T^ ' 

where (j)B{k+) is the leading twist B meson light-cone wave function, normalized according 
to /g°^ dA;+0B(^+) = 1- For the p meson we kept only the twist-2 chiral-odd wavefunction 
0_l('u) 1 16] appropriate for a transversely polarized meson. Since we consider all form factors 



to be renormalized at yU = mb, we will be using this scale in the numerical evaluation 



of (H^). Using fe = 180 MeV, fHrrib) = 144 MeV together with the asymptotic 



expression (f)±{u) = 6u{l — u) and the inverse moment of the B wave function ((A;+) ^) = 3 
GeV-\ we find V{Ep) ~ -(0.9 GeY)/Ep. 

We turn now to a brief discussion of the implications of these relations. The most impor- 
tant result following from (|10|), (|Tl|) was the possibility of extracting the radiative form factor 



g+iy) from semileptonic decays data. Such analyses have been presented in ||T0|, O, O. This 



relation can be now extended to next-to-leading order in A/rnt by combining (pOD with (|TT]) 

g+{y) = -{niB - K)g{y) - V{y) + 0{ml^'') (48) 



Note that including the subleading term in ( pOD is essential, since it is formally of the same 
order in A/ruf, as the other terms kept. 

Numerically, the subleading terms in ( pOD can be significant. For illustration, consider 
the D K* form factors at ?/ = 1, for which numerical results are available from the 
E791 Collaboration ^(1) = -(0.49 ± 0.04) GeV-\ /(I) = (1.9 ± 0.1) GeV fT^. From (gSp 
one finds g^-^^'il) = (0.74 ± 0.06), where we used A = 350 MeV and V{1) = 0. The 
subleading correction in (pO|) contributed about 25% to this result. On the other hand, 
the perturbative effect of hard gluons enters only at two-loop order (through ni, given in 
Eq. (^3D), and therefore can be safely neglected. 

Assuming that the vector form factor g^^^{y) has been measured, the relation (48) 
allows the extraction of g^~^^{y) with accuracy of order A^/m^. Neglecting experimental 



c ■ 

errors, the theoretical uncertainty in such a determination of g^^^{y) could be made as 
small as 10%. A precise knowledge of this form factor is relevant for a good control of the 
short- distance matrix element in rare D decays ^9 . 
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The corrected symmetry relations we derived in this Letter can be expected to help 
reduce the theoretical uncertainties in constrained extractions of form factors from exper- 

1^, and in determinations of CKM 
We will present 
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imental data, such as those presented in []T0|, |Tl] 
parameters with the help of SU(3) and heavy quark symmetry 
details of such an investigation in a separate publication. 
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